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SECOND-ORDER EDGE-PENALIZATION IN 
THE AMBROSIO-TORTORELLI EUNCTIONAL 

MARTIN BURGERi, TERESA ESPOSITO^, AND CATERINA IDA ZEPPIERP 


Abstract. We propose and study two variants of the Ambrosio-Tortorelli functional where the first- 
order penalization of the edge variable v is replaced by a second-order term depending on the Hessian 
or on the Laplacian of u, respectively. We show that both the variants as above provide an elliptic 
approximation of the Mumford-Shah functional in the sense of F-convergence. 

In particular the variant with the Laplacian penalization can be implemented without any difficulties 
compared to the standard Ambrosio-Tortorelli functional. The computational results indicate several 
advantages however. First of all, the diffuse approximation of the edge contours appears smoother and 
clearer for the minimizers of the second-order functional. Moreover, the convergence of alternating 
minimization algorithms seems improved for the new functional. We also illustrate the findings with 
several computational results. 

Keywords: Ambrosio-Tortorelli approximation, free-discontinuity problems, F-convergence, variational 
image segmentation. 
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1. Introduction 

The minimization of the Mumford-Shah functional [31] is one of the most prominent and successful 
approaches in mathematical image processing, which has received considerable attention from a practical 
(see e.g. [5, 11, 34]), computational (see e.g. [20, 25, 34]), as well as theoretical point of view (see e.g. 
[26, 30] and references therein). This approach led moreover to several fruitful variants of the original 
model such as the Chan-Vese functional [9] or region-based variants with realistic noise modelling [33]. 

According to the Mumford-Shah approach and to its later development proposed by De Giorgi and 
Ambrosio [15], the relevant contours of the objects in a picture are interpreted as the discontinuity set 
Su of a function u approximating a given image datum and minimizing the functional 

MS{u) = a [ \yu\'^dx + pn^-\Su), (1.1) 

Jq. 

among all functions u belonging to the space of special functions of bounded variation SBV{Cl). (The 
constants a,/3 > 0 in (1.1) are tuning parameters.) The minimization of (1.1) leads to a so-called free- 
discontinuity problem which is notoriously difficult to be solved numerically in a robust and efficient way. 
An alternative approach frequently used to bypass these difficulties is to relax the minimization problem 
as above replacing the Mumford-Shah functional by a sequence of elliptic functionals which approximate 
MS in a suitable variational sense. A first approximation of (1.1) was studied by Ambrosio and Tortorelli 
[3, 4] who considered a sequence of functionals reminiscent of the Cahn-Hilliard approximation of the 
perimeter. More precisely, in [4] the authors introduced the family of elliptic functionals 

ATs{u,v) = a J {v'^ + r]s)\Vdx + ^ j —helVn]^^ dx, (1.2) 
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with £ > 0 and 0 < ? 7 e <C £, defined on pair of functions u,v G uVit G L^(r2;]R"), where the 

additional variable v encodes an approximation of the discontinuity set of u. Indeed, loosely speaking, 
one expects that a minimizing Vg is close to 1 in large regions of and it deviates from 1, being close 
to 0, only in an £-neighbourhood of Su (where Vug tends to be very large). In this way we get that 
Vg approaches I — xSa as e —>■ 0. This heuristic argument is made rigorous in [3, 4] using the language 
of F-convergence and following earlier ideas developed by Modica and Mortola in the seminal papers 
[28, 29], 

Clearly the above approximation of the Mumford-Shah functional is not the only possible one and 
in particular different variants of (1.2) can be considered in order to enhance the Ambrosio-Tortorelli 
approximation scheme. In this perspective, a computational and practical improvement to the existing 
scheme would be desirable in some cases, e.g. when it is difficult to compute global minimizers of the 
original functional. These considerations also motivate the analysis carried out in the present paper. 
More precisely we are interested in replacing the term e|Vup in (1.2) by a second-order term depending 
on the Hessian or on the Laplacian of v. These second-order penalizations are strongly related to some 
second-order Cahn-Hilliard-type functionals used to approximate the perimeter [17, 24] (see also the more 
recent [10, 12]). At a first glance a higher-order approximation seems counterintuitive, since one expects 
convergence to the first-order perimeter term in the limit. However, we shall observe in computational 
experiments that the second-order approximation may have several advantages. First of all the stronger 
smoothing behaviour of the second-order term leads to smoother approximations of v, which can lead to 
increased robustness in practice. In particular certain structures that are larger than typical noise but 
not yet of interest for the segmentation can be suppressed, e.g. freckles in the segmentation of the face in 
a portrait as illustrated in Figure 1. Moreover an increase in robustness is visible in the computation as 
well, standard alternating minimization algorithms appear to converge faster for the second-order model 
and do not get stuck in undesired local minima, which is the case for the first-order model in a significant 
number of cases. Another interesting aspect is that due to the missing maximum principle in higher 
order equations the optimal value of the variable v is not bounded between zero and one anymore. In 
particular locations where v is larger than one can certainly be identified as edges and due the specific 
shape of the optimal profile one can build two-sided approximations of the edge set in several cases. 



Figure 1. Sisse image, courtesy of Spren Udby (from left to right): Image g, resulting 
V in the Ambrosio-Tortorelli model, resulting v in the second-order model. 

Then, assuming now that v G IF^’^(H) (and dropping the “lower-order terms”) we consider the 
functionals 

Fgiu.v) = a J v^\\7u\^ dx + J —|-£^|V^z;p^ dx, 

as well as the functionals 

E^{u,v) = a J v‘^\Vu\‘^ dx + J —|-£^|Aup^ dx, 
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for which we additionally assume (n — 1) € We show that when e —>• 0 both the families as 

above approximate MS in the sense of F-convergence (see Theorem 4.2 and Theorem 5.3 for the precise 
statements). 

In both cases the most delicate part in the proof is the lower bound inequality. The one-dimensional 
case is first considered and as for the Ambrosio-Tortorelli functional it turns out to contain the main 
features of the problem. In dimension one and Eg clearly coincide (if we ignore the boundary condition, 
which in this case plays no role in the lower bound). When we estimate from below the energy contribution 
of a sequence {ug,Vg) with equi-bounded energy Fg (or Eg), we first appeal to the Gagliardo-Nirenberg 
interpolation inequality to obtain the necessary a priori bound on the first derivative w'. In its turn, 
this bound allows us to deduce that the limit m is a piecewise Sobolev function, which thus has a finite 
number of discontinuities. Then the main difference with respect to the Ambrosio-Tortorelli analysis is 
that now the so-called “Modica-Mortola trick” cannot be applied. Therefore to find the minimal cost of 
a transition between 0 and 1, and occurring in an e-layer around each discontinuity point of u, a careful 
analysis is needed (see Theorem 3.1 and Theorem 4.1). Then, the upper bound inequality follows by an 
explicit construction which yields a recovery sequence (ug,Vg) satisfying (ue — 1) € Wg’ . 

In dimension n > 1 the two functionals Fg and Eg are in general different and the lower bound 
inequality is established in the two cases by means of two different arguments. Specifically, the form of 
the functionals Fg (and in particular the presence of the full Hessian) makes it possible to use a well-known 
integral-geometric argument, the so-called slicing procedure, which allows us to reduce the n-dimensional 
problem to the one-dimensional case and hence to conclude. On the other hand, when we deal with Eg we 
cannot exploit the one-dimensional reduction argument as above because of a symmetry breaking due to 
the presence of the Laplacian. Therefore in this case a different procedure based on the blow-up method 
of Fonseca and Muller [18] is employed. In Theorem 5.3 (see also Proposition 5.1) we show, however, that 
the symmetry breaking at e > 0 disappears in the limit. This is done again using Gagliardo-Nirenberg 
interpolation and then appealing to standard elliptic regularity (where the boundary condition on v plays 
a role). Then in both cases the upper bound inequality follows by a standard density argument and by 
explicit construction. Finally, the F-convergence results are complemented with the corresponding results 
about the convergence of the associated minimization problems (Theorem 4.5 and Theorem 5.5). 

2. Notation and preliminaries 

In this section we set a few notation and recall some preliminary results we employ in what follows. 

Let n > 1; if not otherwise specified, throughout the paper fl C H." denotes an open bounded set 
with Lipschitz boundary. The Lebesgue measure and the fc-dimensional Hausdorff measure on R" are 
denoted by £" and 7^^, respectively. The scalar product oi x,y G R” is denoted by {x,y) and the 
euclidean norm by \x\, whereas A ■ B denotes the product between two suitable matrices A,B. For each 
u G := {x G R”: |a;| = 1}, Q‘' denotes the open unit cube centered at the origin with one face 

orthogonal to u; if xg G R" and p > 0, then Q^(xo) '■= xg + qQ''■ If u belongs to the canonical basis of 
R", we omit the dependence on v and we simply write QQ{xg) = Xg + qQ, with Q := (— 5 , ^)”- 

Let Adh(O) be the set of all bounded Radon measures on H; if p G A4b(f2), we say that pk 
weakly* in A4b(0) if 

/ ipdpk^ / if dp for every (p G Cq (ft). 

Jn in 

Let 1 < p < -too and k G N, we use standard notation for the Lebesgue and Sobolev spaces LP{fl) and 
W^’P{ft). For the general theory of special functions of bounded variation we refer the reader to the 
monograph [2]; here we only recall some useful notation and definitions. For every u G SBV(f2), Vu 
denotes the approximate gradient of u, Su the approximate discontinuity set of u, the generalized 
normal to Su, and it’*' and u~ are the traces of u on Su- We also consider the larger space of the 
generalized special functions of bounded variation on ft, GSBV(fl), which is made of all the functions 
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u G L^(ri) whose truncation := (u A M) V {—M) belongs to SBV{il) for every M G N. We also 
consider the spaces 

SBV^iyi) = {uG SBV{n)-. Vm G L'^{n) and n"-^{Su) < + 00 } 

and 

GSBV^{n) = {uG GSBV{il): Vu G L^{n) and n'^-^{Su) < + 00 }. 

We have 

SBV^{n) n L^{n) = GSBV^{n) n L°^{n). 

Since we heavily use it in what follows, we recall here the Gagliardo-Nirenberg interpolation inequality 
(see e.g. [1, Theorems 4.14 and 4.15]). 

Proposition 2.1. Let U be a bounded open subset o/R” with Lipschitz boundary and let Eq > 0. Then 
there exists a eonstant co(£o, U) > 0 such that 

cqE [ \yvf dx < - f v'^ dx + 

Ju £ Ju Ju 

for every £ G (0,eo] and for every v G 

Moreover, we also recall two a priori estimates for the Laplace operator (see [16, Theorem 1, Section 
6.3.1] and [22, Theorem 3.1.2.1 and Remark 3.1.2.2]) that we use in Section 5. 


Proposition 2.2. Let U be a bounded open subset o/K”. Then 

(i) for each open subset V CC U there exists a constant c{U, R) > 0 such that 

||T||vy 2 , 2 (y) < c{U,V) (||Au|]l 2 ([/) + |]u|]l 2 (£/)) , 

for all V G W^’^(C/); 

(ii) if in addition U has G^ boundary, then there exists a constant c{U) > 0 such that 

Ikllw2.2(;7) < c(17)||A?;|]i2([/), (2.1) 

for all V G H Wg’^{U). 


Throughout the paper the parameter e: varies in a strictly decreasing sequence of positive real numbers 
converging to zero. 

Let a, /3 > 0; we consider the functionals Be, £e : L^{Ll) x L^{Ll) —>■ [0, + 00 ] defined as 


Be{u,v) 



iv-l? 




dx if {u, v) G x 

and uVm G L^(fl;K”), 


+00 


and 


a / dx + 


2^2 . 


{v-iy 


otherwise, 


( 2 . 2 ) 


£eiu,v) := < 


+ e^l Aup ) dx if ( m , v) G x 

and vVu G L^(f2; R”), 


+00 


otherwise. 


(2.3) 










SECOND-ORDER AMBROSIO-TORTORELLI FUNCTIONAL 


5 


We also consider the Ambrosio-Tortorelli functionals ATe- x —>• [0,+oo] 

/3 f Hv-ir 


a J v‘^\Vu\'^ dx + ^ j —|-e|Vtip^ dx if (w, n) G x W^’^(r2) 


ATe{u,v) := < 


+00 


and nVu G R"), 


otherwise, 


(2.4) 


and the Mumford-Shah functional A4S: L^{n) x —> [0, +oo] 

a[ \Vu\^dx + f3n^-\Su) iiueGSBV^in) andv=la.e., 
in 


MS(u, v) := 


(2.5) 


+00 


otherwise. 


3. Optimal profile problem 
In this section we study the optimal profile problem 

f /’+°° 

m ■= inf I {if - 1)^ + if")^) dt: f € (0, +oo), 

/(O) = /'(O) = 0, fit) = 1 if t > M for some M > o|. (3.1) 

The constant m represents the minimal cost, in terms of the unsealed, one-dimensional Modica-Mortola 
contribution in (2.2) and (2.3), for a transition from the value 0 to the value 1 on the positive real 
half-line. 

Theorem 3.1. Let m be as in (3.1); then 

m= min|^^“ ((/-I)" + (/")") dU f & <f(0,+oo), /(O) = /'(O) = 0, ^^^Jif) = l| , (3.2) 

moreover, m = v^- 
Proof. Let 

m:=inf|^^“((/-l)2 + (/")2) dt: / G <^(0,+oo), /(O) = /'(O) = 0, Jm^/(t) = l|. 

By solving the associated Euler-Lagrange equation it is easily shown that rh is attained at 

fit) = 1 + V2e~^ cos(^-^ + ^^ (3.3) 

(see Figure 2). Moreover a direct computation yields th = \pi. 

We clearly have fh < m, then to achieve (3.2) it only remains to show the opposite inequality; i.e., 
fn> m. To this end, we suitably modify / so to obtain a test function for m. Let Xi —>■ -foo as i —>■ -foo; 
it is easy to check that 

lim fixi) = 1, lim ffxi) = 0, 

2 ^ + 00 2 —>- + 00 

with / is as in (3.3). 

We introduce the auxiliary function Q : [0, -foo) given by 

giw,z) := inf |y ((s “ 1)^ + (ff")^) dt: g G ^^([0, 1]), 5(0) = w, g(l) = 1, ^'(0) = 2, gfl) = 0 


(3.4) 
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testing Q with a third-degree polynomial satisfying the boundary conditions, one can easily show that 

liin Cl(ri;,z)=0. (3.5) 

(u;,z)^.(l,0) 

Fix 77 > 0 and let 5 be a test function for Q{f{xi), f'{xi)) such that 

[ {{9 - 1 )^ + ( 5 ")^) dt < G{f{xi),f{xi)) + fj. 

Jo 

Set gi(t) := g{t — Xi) and define 


{ f ( t ) if 0 < t < Xi , 
gi { t ) if Xi < t<Xi + l , 

1 if t > -I- 1. 


Then fi is admissible for m and we have 

r*+oo 


m = 


C “t~00 pXx 

JO JO 

^+00 rXi-\-l 

/ {{h-lf + wmdt- {{9^-l? + {9'lf) 

JO Jxi 


dt 


> m - g { f { xi ), f '{ xi ))- r ]. 

Invoking (3.4) and (3.5), we conclude by first letting i —>• -l-oo and then 77 —?> 0+. 


□ 



Figure 2. The optimal profile /. 


Remark 3.2. For d G R, set 

nid := minj^ ((/- 1)^ -k ( ff ) dt : f G (0,-koo), /(O) = d , /'(O) = 0, ^Jfm^/(0 = l| • 
A direct computation gives 


hence lim(;_>o rrid = y/2. 


nid = V2{d- 1)^, 


(3.6) 
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4. The first model 

In this section we study the F-convergence of the functionals defined in (2.2). The one-dimensional 
case is considered first: this is the object of Subsection 4.1. Then, Subsection 4.2 is devoted to the 
n-dimensional case, with n > 2. Finally, in Subsection 4.3 we deal with the equicoercivity of a suitable 
modification of and with the convergence of the associated minimization problems. 

In all that follows c denotes a generic positive constant which may vary from line to line within the 
same formula. 


4.1. The one-dimensional case. Let —oo < a < b < -|-oo and let L^{a,b) x L^{a,b) — >■ [0,-|-oo] 
be the one-dimensional version of the functional in (2.2); i.e., 


^e{u,v) 


-foo 


if {u, v) G IT^’^(a, b) x lF^’^(a, b), 
otherwise. 


Notice that in this case the condition vu' G L'^ia.b) is automatically satisfied for any pair (u.v) G 

W^-^{a,b) X lT2-2(a,6). 

We have the following F-convergence result. 


Theorem 4.1. The sequence (J>) T-converges, with respect to the {L^{a,b) x (a, b)) -topology, to the 
functional AiS : L^{a,b) x L^(a,b) —> [0,-foo] given by 


A4S(u, v) 



J a 


-foo 


ifu G SBV^{a, b) and v = 1 a.e., 
otherwise, 


where ff{Su) denotes the number of discontinuity points of u in (a,b). 


Proof. For the sake of notation, for any open set U C (a, b) we consider the localized functionals 


Ae{u,V, U) 



if {u,v) G X IT2'2([7), 

otherwise. 


We divide the proof into two steps in which we analyze separately the liminf-inequality and the limsup- 
inequality. 

Step 1: liminf-inequality. Let {u, v) G L^(a, b) x L^{a, b) and (ug, Ug) C L^(a, b) x L^(a, b) be such that 
{ue,Ve) {u,v) in L^{a,b) x L^{a,b). We want to prove that 


liminf/g(ug,Ug) > M.S(u,v). 


(4.2) 


Clearly it is enough to consider the case 

liminf/g(ug, Ug) < -foo. 
£^0 


Then, up to subsequences, Ijug — l\\L^(a,b) < c^/e, from which we immediately deduce that v = 1 a.e. in 
(a, 6). 

According to Proposition 2.1 there exists a positive constant cq > 0 such that for e sufficiently small 
we have 

CO feiv',)^dt< f'e^v'ffdt. 

J a J a ^ J a 
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Therefore 



dt < c 


for e > 0 sufficiently small; 


hence we can apply [7, Lemma 6.2 and Remark 6.3], with Z = {1} and W{s) = (s — 1)^ to conclude that 
there exists a finite set S such that, for every fixed open set 17 CC (a,b) \ S and for e: > 0 sufficiently 
small, 1/2 < Ue < 3/2 on U. For every such fixed U we have 


^ sup / {u'^Y dt < a sup / dt < +oo; 

4 £>0 Ju £>0 Ja 


thus u G W^’'^{U) and Mg ^ m in Moreover, since Ug —>■ 1 in L^(17) and u'^ u' in L'^{U), then 

u'^Vg u' in L^{U); hence we have 


a f (m')^ dt < a liminf f dt < a liminf [ vtiu'S^dt. 

Ju ^^0 Ju e^o 

By the arbitrariness of U, (4.3) can be rewritten as 


(4.3) 


' (a,b)\{S+l-rj,rj]) 


{u'Y dt < 


dt < a liminf 

e ->0 


(a,f))\(S+[-77,77]) 


:{<fdt< 


dt < a liminf 
£— ^0 


f 


Ku'e^dt, 


(4.4) 


for every r] > 0. This allows us to conclude that u G SBV‘^{a, h) and Su C S. 

Let N ■= Su ■= {ti,... ,tAr} with ti < t 2 < • • • < tjv, and consider pairwise disjoint intervals 

li = (tti, bi) C (a, b) with U G li] we want to show that, for every i G {1,..., N}, 


liminf J'e(M£,Ue,Ii) >/3. (4.5) 

£->■0 


To this end, fix t G {1,...,A^}. Let ti G /' CC h and set mt := liminf£_>o inftg// (z;E(t))^. Since /' 
contains a discontinuity point of u, we want to show that it is convenient for to be “close” to zero in 
i.e., that = 0. We argue by contradiction assuming that > 0. Then for every t G we have 


lim inf 
£— ^0 


nil 


> 1 , 


therefore 


[ {u'^fdt<—f vl{u'^fdt< 

J I' nii j j! 




TUi 


< c. 


SO that Me ^ u in VF^’^(/') and Su H /' = 0 which contradicts the hypothesis U G Hence, we must 
have nii = 0. As a consequence there exists a sequence (s®) C /' such that 


lim Me(4) = 0 . 

e —>-0 


(4.6) 


On the other hand, up to subsequences (not relabelled). Me —>■ 1 a.e. in (a, 6) as e —>• 0. Therefore, there 
exist f®, r® G A such that f® < s® < r® and 

lim Me(r®) = 1, lim Ve{r\) = 1. 

£—>•0 £—)-0 


(4.7) 
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Moreover, again appealing to the interpolation inequality Proposition 2.1 we may deduce that eu' —^ 0 
in L^(a, h). Indeed 

\ 1/2 


J e\v'^\dt < e (u' 


< e" 


+ e^v':y dt 


1/2 


< —>■ 0 as e —>■ 0. 

Then up to subsequences (not relabelled) it is not restrictive to suppose that 

lim£u'(f*) = 0 and lim eu'(r*) = 0. 


E —^0 


E —^0 


Let s* be a minimum point for in In view of (4.6) and (4.7) we deduce that s* G for 

e sufficiently small. Then, since £ C^([f*,r£]) we have v'^{sl) = 0. Moreover Ue(s*) < Vs{sl) < ^^(r*), 
hence (4.7) implies that Ve{sl) < 1 for e small. 

We have 

(''-^^^+£3(^;")2^ dt 


(He I’^ei ^i) — 


/3 


> 


2^2 

P 

2^2 




We now estimate from below the term 


Jsi lo ~ 

where z = {t — s\)/e and We{z) := Ve{ez + s®). To this end let ge,i £ C^([0,1]) be an admissible function 
for g{veirl),ev'^{rl)); i.e., g^diO) = Ue(r*), 5 e,i(l) = 1, g',(0) = sv'Jri), g'_i(l) = 0. By construction 


lim g^ i(0) = 1 and lim g' hO) = 0, 

£->0 ’ ES-O 


hence by (3.5) we infer 


lini^;(ue(r*),£u'(r*)) = 0. 


£->0 


(4.9) 


Let (ue,i) be the sequence defined as 

We(z) 

e,i ( 2 ^) ■— ' ge,i ^Z ■ 
1 


if 0 < z < 


K - K 


if ^^ < z < ^^ + 1, 


if z > ^ + 1. 


By definition of ge,i it follows that {ve,i) C Wjq’^(0,+ 00 ). Since (u£,i) is a test function for (with 

as in (3.6) with d = Veisp))^ we have 

r+^ 

/ {{We - if + {w'^f) dz = / {{Ve,i-lf + {Ve,if) dz-g{ve{rl),evfrl)) 

Jo Jo 

> m^,{si)-Q{ve(rl),ev'^{rl)). 
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A similar argument applies to the term 


i: 

J rl 


(Ve - 1 )^ 


+e\v':r dt 


in (4.8). Since = 'J2{ve{s\) — 1)^ and Ve{s\) < Ue(s*) < 1, we have 

Te{ue,Ve,Ii) > - ^g{ve{rl), E v'^{rl))) > ^^inin. ^ £ <(r*))) 

^ = ^{v,{sl) - if - -^g{vfrl),ev'frl))) 

and (4.5) follows from (4.9) letting e —?► 0. 

Since the intervals A are pairwise disjoint, gathering (4.4) and (4.5), we get 


/3 

liminf J>(m£, lie) > a liminf j v^{uf'^dt-\ -^ liminf 


£—^0 


£—^0 


2^2 


iminf / 

"^0 Ja V 


{Ve - 1 )^ 


+e%v':r dt 


> a 


iu'fdt + P^iSu). 


J {a,b)\{S+[-ri,ri]) 

Hence the liminf-inequality (4.2) follows from the arbitrariness of rj. 

Step 2: limsup-inequality. Let u G SBV^{a, 5), Su := {ti,..., t^v} with ti < t 2 < ■ ■ ■ < tN, and set 
to := a, tAT+i := 6, Sq ■= minj{tj_|_i — : i = 0,..., N}. For i = 1,..., iV, define A := , *‘~*~*‘~*'^ . 

We now construct a recovery sequence for the F-limit. 

Let r] > 0; there exist a function G Wjof (0, +oo) and a constant Mn > 0 such that fniO) = /))(0) = 0, 
fri{t) = 1 for all t > Mri, and 

rM„ 

{{fr,-if + {f:;f)dt<V2+7j. (4.10) 

Let > 0 be such that —>■ 0 as e —>■ 0, then for e sufficiently small For i = 1,..., iV let 

ifl be a cut-off function between (ti — and {U — + ^e); *.e., ipl G C^(ti — + ^e), 

0 <</?£< 1 and = 1 on (ti — Define the sequence 


/ 

^0 


fx) := u{x) ^ ^''eix)xii (a;)^ ; 


then [ue) C h) and Mg —>■ m in L^(a, h) by the dominated convergence Theorem. 

Fix T > Mr)] then > T, for e small. Define the sequence 


fO if \t-ti\ <^£, 

'*^e(^) := I if < |f - L| < ■?£ +£T, 

if if t€{{\t-tr\>^e + eT}nlr)u[a,^]u[i^,b], 


for f = 1, • • • , A^, where fr) is as in (4.10). It is immediate to check that (ug) C VF^’^(a, b) and Ug —>■ 1 in 
L^(a, b). Moreover for every i = 1,..., A^, we have 



L 


Ii\[ti-is,ti+ie] 


]{u'f dt 


(4.11) 
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and 


L 

L 


di 


ti-ie ( 1 

-\u 






- 1 ) 




dt 


■L 


'ti+^e+^T' 


ti+i 

T 




3 I ft! 
V 


= 2 


[ iifvi^) - 1)^ + (4'(z))^) dz + 
Jo 


- 1 ) +e^f, 

e 




dt ■ 


rU+is 1 

4-€e £ 


dt 


< 2 V 2 + 7 J + 2 — 


(4.12) 


Therefore gathering (4.11) and (4.12) gives 


d' ^e) — 


N N 

^ I vUuTdt+-y=^ 




N 


^ «E/ 

i=lJli\lti-ie,U+ie\ 

Finally, invoking the dominated convergence Theorem yields 


vliu)'^ dt + f}N + c— + 77 . 

e 


pb 

limsup J'E(Me,i;e) < a / (u')^ dt +/3 #(S'„) + 77 

£—>^0 J a 


and thus the thesis. 


□ 


4.2. The n-dimensional case. Let n> 2 and let be the functional defined as in (2.2). The following 
F-convergence result holds true. 

Theorem 4.2. The sequence (d>) T-converges, with respect to the x L^{fl))-topology, to the 

functional A4S: x L^(n) —> [ 0 ,+oo] given by 

, , ia f \Vu^\dx +^n'^-^(Su) ifueGSBV’^iil) andv = l a.e., 

7W5(u,?;):=< Jn 

[ +00 otherwise, 

Proof. As for the one-dimensional case, if A is an open, A C fl, we set 

/3 f f{v-ir 


aj v‘^\Vu\'^ dx + ^ J + e3|y2^|2^ A {u,v) €W^-^{A) xW^'^{A) 


Be{u,v,A) := < 


-|-oo 


and 


MS{u, V, A) := < 


and vVu € L'^{A-W^), 
otherwise, 

af \\7u\^dx + l3n'^-\Sur\A) iiu £ GSBV^{A) and v = 1 a.e., 


^ +00 otherwise. 

We divide the proof into two steps in which we analyze separately the liminf-inequality and the limsup- 
inequality. 
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Step 1: liminf-inequality. Let A be an open subset of Q. We recover the lower bound 

F- liminf A) > j\4S(u, 1, A) (4-13) 

e ->0 

from the one-dimensional case, by using the slicing method. To this end we start recalling some useful 
notation. For each ^ G we consider the hyperplane through the origin and orthogonal to 

n« := {x e R”: (x ,0 = 0 }, 

and, for every y € 11^, we consider the one-dimensional set 

A^^y != {t € R: € A^ 

and the one-dimensional functions defined by 

u^,y{t) := u{y + t^), v^^y{t) := v{y + t^). 

By Fubini’s Theorem we have 


Te{u,v,A) 

= f f {av^{y + t^)\Vu{y + t^)\^ 

Jni J A(:^y \ 

> f f (av^{y + t^)\{Vu{y + t^),0\^ 

Jni JAf^y V 

^ f ~ + e^iivMy + tot op) ) dtdn^-\y) 

Jn( 


where for all ^ and y, are the one-dimensional functionals defined as in (4.1). 

Let (m,u) G L^(f2) X and let C L^{il) x L^(f2) be such that {u^,v^) -A {u,v) in 

L^{n) X L^(n), and 

liminf J>(u,i;, A) < -boo, (4-14) 

£—>•0 

then, V = 1 a.e. in A. By Fubini’s Theorem and Fatou’s Lemma, 


0 = lim / Irte — u| dx 


im [ 

lim / [ IM^^y-u^Jdtdn^-^y) 
7n« Ja, y 


> 


/jjj Il(«e)«,y - H,y\\L^{Ag.y) d'H'' (y); 


hence —>■ in L^{A^^y) for "H" ^-a.e. y G and analogously {ve)^,y -A 1 in L^{A^^y) for 

7^"“^-a.e. y G 11^. Therefore, appealing to the one-dimensional result Theorem 4.1 we have that G 
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> 


> 


SBV‘^{A^^y), for "H" ^-a.e. ?/ G 11 ^ and 

\iminiB^{ue,Ve,A) > liminf [ {{ue)^,y, ive)^,y, A^^y) 

£->■0 e - s-0 Jyis 

[ liminf y, A^^y) 

a\uly\^dt + p#iSu,,,nA^,y)^ d-H--\y). (4.15) 

Let M G N and consider the truncated functions um '■= {—M) V (m A M); by (4.14) and (4.15) we have 

f \DiuMh,y\iA^.y)d-H--\y) 

Jn( 

= l('«M)^,y| dt + |('«M)^y - (MM)^y|#(5'(„M)j,„)^ '^^"”^(2/) 

< \Mlyfd?j +2M#(5(„,)^ jj d7^"-i(y)<+oo. 

Thus, by virtue of [ 8 , Theorem 4.1(b)] we conclude that um G SBV'^{A) for every M G N, hence the 
lower semicontinuity of the Mumford-Shah functional entails u G GSBV'^{A). Moreover, by (4.15) and 
[8, Theorem 4.1(a)], we infer 

L-liminf 7e(u, w, A) > a f \{Vu,^)\‘^ dx + P f j(^, z/„)j (4-16) 

'^“*'0 Ja J AnSu 

In particular, since (4.16) holds for each ^ G S"“^, we have 

T-liminf J^e(M, u, A) > sup la f \{Vu,^)\'^ dx + f3 [ 1 (^, u^)] 1 . 

JGS"-1 I JA JAnSu J 

The previous estimate can be improved to (4.13) by means of a measure theory lemma. Observe that for 
u G GSBV‘^{A) and v = 1 the set function y.{A) := L-liminf£_>o .Fe(m, n, A) is superadditive on disjoint 
open sets. Then, applying [7, Lemma 15.2] with 


A:= £”+'H”"V5„ 


and 


al(VM,^i)l^ ifx^S'^, 
P\{^i,i'u)\ iixGSu, 


ipi{x) ■■= 

where (^i) is a dense sequence in §"“^, we obtain that 

T-liminf J>(m, n, A) > / sup'^idA 

^“*■0 J A i 

= ctj \yu\^dx + pn^-^{Sur\A) 

= A4S{u,v,A). 

Finally the liminf-inequality follows taking A = fl. 

Step 2: limsup-inequality. We want to show that 

T- limsup J>(u, v) < A4S(u, v), 

6 —^0 


(4.17) 
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whenever u € GSBV^{^1) and v = 1 a.e. in il. To this end it is useful to recall the density result [13, 
Theorem 3.9 and Corollary 3.11]. Let u £ GSBV^{il) then there exists a sequence (uj) C SBV'^{il) 
satisfying the following properties: 

( 1 ) Suj is essentially closed; i.e., n \ Suj)) = 0 ; 

( 2 ) Suj is the intersection of with a finite number of pairwise disjoint closed and convex sets, each 
contained in a (n — l)-dimensional hyperplane, and whose (relative) boundaries are C°°; 

(3) Uj £ \ Suj) for every fc € N; 

(4) \\uj - m||l 2 (q) -)> 0 as j -)■ +oo; 

(5) \\Wuj - Vm||l 2 (q) 0 as j +oo; 

( 6 ) \n^-\Su,) - n^-\Su)\ ^ 0 as J ^ + 00 . 

Denote by W(D) the class of all functions for which conditions (l)-(3) hold. Then, by a standard density 
argument it is enough to prove (4.17) when u belongs to >V(D). Indeed assume (4.17) holds true in >V(D). 
If u £ GSBV‘^{V,) then there exists a sequence (uj) C >V(D) satisfying (4)-(6). Hence it holds 


L-limsup J>(it,i;) < liminf 

£->0 t-S'+oo 



lim sup Jg {uj , 
6 —^0 



< lim A4S(uj,v) = A4S(u,v), 

j—^-hoo 


where the first inequality is due to the lower semicontinuity of the L-limsup, whereas the last equality 
follows from (5) and (6). 

We now prove (4.17) for a function u G W(r2). Assume first that Su = fid K where AT is a closed 
convex set contained in an (n — l)-dimensional hyperplane H with normal v. Let p: R" —>■ H be the 
orthogonal projection on H, d{x) := dist(a;,n), and for any (5 > 0 set := {x £ H: dist(a;, AT) < (5}. 

Let ry > 0; then there exist a function £ WjQf(0,+oo) and a constant > 0 such that /,,(0) = 
/'(O) = 0, fnit) = 1 Vt > and 



+ dt<V2 + r,. 


(4.18) 


Fix T > Mjj and let > 0 be such that —>■ 0 as e —0; set 


Ag := {x £ R”: p{x) £ A'®, d{x) < 

Bg := {x £ R": p(x) £ , d{x) < ^g + eT} 


(see Figure 3). Consider moreover 



Gg := {x £ R": p(x) £ A:'^/^ d{x) < Ce/2} 
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and let be a cut-off function between and A^; i.e., G 0 <(/?£< 1 and = 1 on C^- 

Define the sequence 

Ue{x) := M(a;)(l - Pe{x))\ 

then {ue) C and Mg —)► u in L^(D), by the Lebesgue dominated convergence Theorem. 

Let 7 e be a cut-off function between and i.e., G C^{K‘^^), 0 < 7 ^ < 1, = 1 on with 


l|V7slU=^(n) < 

£ 

||VSlU~(n) < 

for some c > 0. Let he: K. —R be the function defined by 

0 if \t\ < ^e, 

1^1 - e 


(4.19) 

(4.20) 



where frj is as in (4.18). Let 


if Ce < \t\ < Ce + e'T, 
if \t\>^,+eT, 


Ve(x) := -ie{p{x))he{d{x)) -h (1 - 7 e(p(x))); 


(4.21) 


by construction (ve) C IT^’^(D) fl L°°{id) so that UeVrte G L^(0;]R"); moreover Ug —>■ 1 in L^(0). 
We have 


limsupJ'e(we,r’e) = limsup / (aUel^^^EpH—^ f— -—he^lV^iiep 


e->0 


£—>^0 J n 


2^/2 


6—^0 J n 




< a limsup / v^lVuef dx -\ - -j= limsup 


2^2 


£—^0 


■/..( 


dx 


+g3|^2^ | 2 \ 


—-= limsup / 

2v2 £—^-0 J Bg\Ag 


I .-3|v72„, |2 


-h£^|V"Uer dx 


lim sup 


In\Be ( 


2-\/2 e^o Jq\b, 




=: + /s + -^ 4 )- 

We now compute separately the four terms Ii, I 2 , 13 , and I^. 

Since Ug = 0 on Ae and Ue = m on D \ Ag, we immediately deduce 


h = 


limsup / v^\Vu\'^ dx = / iVwpda;; 

e-lO Jn\Ae Jo, 


(4.22) 


(4.23) 


moreover the term I 2 gives no contribution to the computation, indeed 


I 2 = lim sup 


-f 

e Ja, 


dx 


e^O £ J A 

1 B 


= lim sup - f dt f d'H^ ^ (y) 

£—^0 ^ J — J J<C^ 

= 21 imsup —■H”-i(iL'^) = 0 , 

£->0 e 


where we have used the fact that >C £ and ^(K^) -A ’H" ^{K) < -l-oo as £ —>■ 0 . 


(4.24) 
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Since = 1 on \ Bg, we also have 1/^ = 0. Then it only remains to compute I^. To this end it is 
convenient to decompose \ as the union of two sets, Dg and iSg, defined as follows 

T>g := {x G M”: p(x) G K^, ^g < d(x) < Ce + sTj 
E, ■= {x G M”: p{x) G \ K^, d{x) < Ce + eT} 


(see Figure 4). 



Figure 4. The sets Dg and E^. 


On £)g the function Ug = /j, moreover V(i(a;) = ±u hence we have 

f 

J D, 


lim sup 
s —>-0 J 


+ e^\V^v. '2 


= lim sup 

€—^0 J Dp 


fri 

^e-\-eT 


d{x) - ^g 


< 2 lim sup 

£->■0 Jk ^ J (, 


LI a 


dx 

- 1 1 + 
t-C 


d{x) - v'^ -V 


- 1 ) +-(/: 


t-ie 


dx 


dtdn^-^iy) 


= 2 1 imsup/' f {{fr,{t)-lf + {f”{t))^)dtd'H'^^{y) 

£->■0 JJo 

< (2^2 + ^) limsup'H”-^(iG") = i2V2 + fi)n^-\K). 

£—^0 


(4.25) 


Hence finally to achieve the limsup-inequality we have to show that 

limsup f f— -—|-e^|V^Ug|^^ dx = 0. (4.26) 

£->■0 Je^ V £ / 

To do this we consider the further decomposition of Eg as the union of two sets 14 and Wg dehned by 

14 := {x G K": p{x) G \ , d{x) < ^g} 

Wg := {a; G K": pix) G \ , Ce < d{x) < ie + eT} 

(see Figure 5). Since he{d{x)) = 0 on 14, we have Ve{x) = 1 — 7 g(p(a;)) for x G 14- Then, if we denote by 
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Figure 5. The sets K and W^. 


Dp(x) the Jacobian matrix of p evaluated at x, we get that ||iJp(x)||ioo(Rn.jjnxn) < 1 and 

(Vg — 1 )^ ^ 3 |v 72 „, |2^ 


i(‘ 


+ e^\V^VeV dx 


IVe 




+ \-Dp{x)^ ■ Dp{x)\ ) dx 


< 


/ / 


^ + £3|vS(2/)r) dtdn^-\y) 


< c — n^-\K^^\K^)^0 ase^O. 
e 


where in the last inequality we have used (4.20). 
Also 

limsup / 

6 —>^0 J 


l!i_Jl+e3|y2^^|2 


dx = 0, 


indeed Ve{x) = Je{pix))fp + (1 ~ leipix))) on Wg and we have 


/ 

JWe 


) dx 


J 


= 2 




Jt 

+e^ Dp{y + tvf^ ■ y'^le{y) ■ Dp{y + tv) ( /, 


- le{y) 

t-i 


- 1 


44^ 


{4tle{y) ■ Dp{y + • Vd{y + tv) + Vd{y + tv)'^ ■ V'-feiy) ■ Dj 


M 4 


t- ^e\ v'^ ■ V 


dt dH^ (y) 

n—1 ( T^2e ' 


< -e'H^-^{K^^\K^) = c'H'^-^{K^^\K^)^Q as £ ^ 0, 

e 

where to establish the last inequality we have used (4.19) and (4.20). 


{y + tv)) 
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Thus gathering (4.22)-(4.26) we finally deduce 

limsup Ue) < a f \\7u\'^ dx + {/3 + fj)(Su) V 77 > 0 , 
s—^0 J Cl 

and hence the limsup inequality. 

We now consider the general case in which Su = with Ki closed and convex set contained 

in an (n — 1)-dimensional hyperplane 11®, with normal Let pi : R" —>■ 11® be the orthogonal projection 
on n®, di{x) := dist(x, 11®), and for every <5 > 0 set K- := {x G If®: dist(x, Ki) < 5}. 

Let frj be as in (4.18) and fix T > for 0 < ^e -C e consider the sets 

Al := {a: € R®®: p,{x) € Kf, di{x) < ^e} 

B® := {a: G R®®: Pi{x) G , di{x) < -f eT}. 

Arguing as above we can construct a recovery sequence (m® , vl) for Fe{u, v, Ki). Then, let eq > 0 be such 
that 

VO < e < Eo, 


T< 


2e 


(4.27) 


where Jq := min{dist(A'i, ATj ): i, j = 1,... ,r, i j}; for e G (0, Eq) set 

r r 

A,-.= \JAI, Be-=[jBl 

and consider the sequences (ue) C W^’^(n) and (ue) C W^’^(n) fl L°°{n) defined by 


and 


Ueix) := 


-{x) 


u® if a; G A®, 
u if a; ^ Ae, 

u® if a; G B®, 
1 if a; ^ Be- 


(4.28) 


(4.29) 


Condition (4.27) ensures that i?® (1 B^ =0 for all i ^ j and consequently that the sequences in (4.28) 
and (4.29) are well-defined. Moreover we have (ue,Ue) (at, 1) in L^(0) x L^(0) and 


limsup J'e(ue,Ue) 


£—^0 


< a limsup 

£—>•0 J Q 


f vl\V 

JQ. 


UeT dx + 


2 V 2 V 


lim sup 
£^0 J Ae 


{Ve - 1 )^ 


-I- e^IV^iIeP 1 dx 


'L 


■ lim sup 

£—^0 J Be\Ae 


1)^ I ^3|v72„, |2 


-l-E |V UeI dx-|-limsup 


e->0 JQ\B, 


I 

ao' 


(!(^^+e3|v2^ ,2 , 


< a f \Vu\‘^dx-\ —^ ( limsup f 

- Jn 2^2 e^o^JAi 


(Ve - 1 )^ 


-I- e^IV^UeP 1 dx 


'L 


■ lim sup 

£—)-0 J Bl\A\ 


1)^ I ^3|v72„, |2 


-l-E |V UeI dx-|-limsup 


£->0 aa\B, 


/ 

ao' 


^)% g3|v2^|2 , 


a[ \Wu\^dx + 13'^n’^-AK,) 

i=i 

a[ \\7u\^dx +pn^-CSu)=MSiu,l)- 
Jn 


□ 
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Remark 4.3. Let / G VLg’^(r2) then it is immediate to show that 

ll^/l|L2(n) = ||V^/|li2(f2;RnXr.). 

Therefore, if m G (i; — 1) G and vVu G L^(r2;R") the functionals Te can be rewritten 

(u-l)2 


as 


J'e(u, v) = 


) = a f 

Jn 


dx + 


J_ 

2y/2 


+ Auf 1 dx. 


Let J-e- X L^(f2) —>■ [0,+oo] denote the functionals 


a / dx + 


2 V 2 


I 

Jn 


{v-iy 


:Fe{u,v) := < 


+ e^|Aur) dx if{u,v-l)GW^^^{n) xWo'\n) 
and uVu G K”), 


^ +00 otherwise, 

then the following theorem is an immediate consequence of Remark 4.3 and Theorem 4.2. 

Theorem 4.4. The sequence (J>) T-converges, with respect to the (L^(n) x L^{Q))-topology, to the 
functional AiS, defined as in (2.5). 

Proof. Since every sequence (ug) with equibounded energy strongly converges to 1 in L^(n), the liminf 
inequality is trivial. Hence, it only remains to prove the limsup inequality. To this end we notice that 
the recovery sequence exhibited in Theorem 4.2 is also a recovery sequence for the functional J>. Indeed 
it is immediate to check that (ve) defined as in (4.29) (see also (4.21)) satisfies {ve — 1) G ITg’^(r2). □ 

4.3. Convergence of minimization problems. Let ? 7 e > 0 be such that rj^/e —0 as e —0. Let 7 > 0 
and g G Lf{Ti) be given and for every {u,v) G IT^’^(H) x IT^’^(H) consider the functionals 


Te{u,v)-\-rie |Vupfia; + 7 / \u - g\ 
Jn Jn 


- grdx. 


(4.30) 


The following result holds true. 

Theorem 4.5. For every fixed e > 0 there exists a minimizing pair {ug,Vg) for the problem 
Mg := inf I J^e(M, v) + 77e / \Vu\'^ dx-\-j f |u — gp dx: (u, u) G IT^’^(H) x IT^’^(H) 

L J Q, J Q 

Moreover, up to subsequences, {ue,Ve) —>■ (u, 1) in x L^(H) where u is a solution to 

M := min •|a45('«, 1) -\-'j J |u — gp dx: m G GSBV^{il)^ ; 
if n = 1 then u G SBV'^{fl). Further, Me —>■ M as e —>■ 0. 


(4.31) 


Proof. For fixed e > 0 the perturbation term rje Jq |Vup makes the functionals in (4.30) coercive with 
respect to the weak (IT^’^(H) x lT^’^(H))-topology. Indeed, let {uk,Vk) C IT^’^(H) x IT^’^(H) be such 
that 


Ae{uk,Vk) + Pe IVUfcpdx + 7/ \uk-g\'^dx 

Jn Jn 

As a consequence we deduce 

||i^feVMfe||i 2 (n.Rn) < c, ||Mfc||ivu 2 (n) < c 


Me as /c —>• + 00 . 
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for some c > 0 independent of fc; moreover, by the interpolation ineqnality Proposition 2.1 we also have 
||'Pfc||vi/ 2 , 2 (Q) < c. Then up to subsequences (not relabelled) 


(4.32) 

Vk^v in 

Therefore 

VkVuk^vVu in L^(n;R”), 

hence by the uniqueness of the weak limit v'S/u = w G i^(f2;M”). Then, the existence of a minimizing 
pair (ue,Ve) easily follows appealing to the weak lower semicontinuity of the norm and to the direct 
methods. 

The requirement that r^e/e—^-Oaser—>-0 ensures that 


VkTuk w 

in 

Uk ^ u 

in 

Vk ^ V 

in 


J^e(M,w)+r 7 e / |Vup(ia: + 7 / \u — g\'^dx MS{u,v)+'j / \u — gY 


dx 


(4.33) 


Jn 


Jn 


with respect to the strong {L^{n) x T^(f2))-topology. This can be easily seen arguing as in Theorems 4.1 
and 4.2 (now taking e) and recalling that T-convergence is stable under continuous perturba¬ 

tions. 

Moreover, if (Mgjiig) C IT^’^(O) x IT^’^(n) is any sequence such that 


sup i J^e{Ue,Ve) + ge \'VUeYdx+ / \Ue - gY dx 
e V Jn Ja 

then the interpolation inequality Proposition 2.1 implies 

sup i AT e{Ue,Ve) + ge f dt + f \Ue - g\ 

^ \ J a J a 


— q\ dt 


< - 1 - 00 , 


< - 1 - 00 , 


where ATe is as in (2.4). Thus [4, Theorem 1.2] immediately yields the equicorecivity of the functionals 
as in (4.30). Finally, by virtue of (4.33), the convergence of the associated minimization problems is 
ensured by the fundamental property of P-convergence. □ 


5. The second model 

In this section we study the asymptotic behaviour of the energies Sg defined in (2.3). Specifically we 
prove that, up to imposing suitable boundary conditions on n, the P-limit of fg is again given by (2.5). 

In what follows il will be an open bounded subset of R" with boundary. 

Let £e be as in (2.3); we have 

£s{u,v)<a f v'^\TuYdx-\ —f (— -^—|-2e^|V^np^ dx, (5.1) 

Jn 2 V2 Jn\£ J 

for all {u,v) G x IF^’^(fl) such that vVu G L^(f2;M"). 

It is convenient to introduce the following notation. Let £',£"'■ L^(£l) x —>■ [0,-|-oo] be the 

functionals defined as 

^'(•c) := r-liminffg (•,•), £"{■,■) := P-limsupfg (•, •); 

E-i-o g^O 

then, by virtue of Theorem 4.2, we get 

£'{u,v) <£"{u,v) <a [ \VuYdx + 2l3n^-\Su) 

Jn 
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for all u £ GSBV^{Vt) and for v = 1 a.e. in 17, hence 

GSBV^{fl) X {n = 1 a.e. in 17} C domS" C domf'. 

We now apply the blow-up argument of Fonseca-Miiller [18] (see also [6]) to obtain the following lower 
bound inequality for the functionals 

Proposition 5.1. For every u £ GSBV^{El), we have 

r-liniinff£(u, 1) > MSlu, 1), 

e->0 

with M.S defined as in (2.5). 

Proof. Assume first that u belongs to SBV^{El). Let {ue,Ve) C L^(17) x L^(17) be such that (uejUg) —>■ 
(u, 1) in L^(17) X L^(17) and sup^ Ue) < +oo. 

For each e > 0 consider the measures 

h-s ■■= + -^J= £” l 17 . 

By hypothesis /ie(17) = is equibounded therefore, up to subsequences (not relabelled), /Xg /x 

where /x is a non-negative finite Radon measure on 17. Using the Radon-Nikodym Theorem we decompose 
/X into the sum of three mutually orthogonal measures 

^i = ^J^aC^+^iJ%'^-KSu+^Xs 


and we claim that 

h-a{xo) > a |V?x(a;o)|^ for £"-a.e. xo € fl (5.2) 

and 

h-jixo) > fi for 'H"'“^-a.e. xq £ Su- (5.3) 

Suppose for a moment that (5.2) and (5.3) hold true, then to conclude it is enough to consider an 
increasing sequence of smooth cut-off functions (</2fc), such that 0 < < 1 and swpj.ipk{x) = 1 on 17, 

and to note that for every k G N 


lim. SJup.Vp) > 


> 


lim inf 

S —^0 




(IVxxe 


fi ({Ve- 1 )^ 


2^2 V e 

/ Fkdfj.> / ^aPk dx+ yLjipk d'}r‘~^ 
Jn in JSu 

a / \Wu\‘^ipkdx + fi / (pkdH 
in is,,. 


-f e^|Ax 


<Pk dx 


jn — 1 


Hence, letting k —> -foo the thesis follows from the monotone convergence Theorem. 

We now prove (5.2) and (5.3). We start proving (5.2). To this end let xq £ 17 be a Lebesgue point for 
/X with respect to £” such that u is approximately differentiable at Iq; be., 


tJ-a{xo) 


y KQei^o)) 
s^oC-{Q,{xo)) 


h'iQejxo)) 

g^O p" 


(5.4) 


and 

lim—^ [ |xx(a;) — ?x(a;o) — (Vxx(a;o),x — xo)! da; = 0. (5-5) 

iQ,(xo) 

The Besicovitch derivation Theorem together with the Calderon-Zygmund Theorem ensures that (5.4)- 
(5.5) hold true for a.e. xq £ 17. 
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Since ^ is a finite Radon measure we have that ^{dQ^{ xq)) = 0 for all g > 0 except for a countable 
set. Thus, being Xq~(^ upper semicontinuous function with compact support in (for g small) [2, 
Proposition 1.62(a)] yields 


If 1 

fia{xo) = lim — /_d//> limlimsup —^£(Qj,(xo)) 

Q jQJxn) £->0 Q 

av^{x)\Vue{x)\ 


1 f 

= lim lim sup — / 
g_^o P" Jc 


Qe(xo) 


2 , P ({Ve{x) - 1 ) , ^^^,2 


2^2 


+ e'^|Ane(x)| dx 


= limlimsup/ [ avl{xo + Qy)\\/Ue{xo + Qy)Y 
e->o Jq V 


/3 f {ve{xo + gy) - 1) , ^3 


2^2 


+ e:'^|Ang(xo + gy)\ dy. 


Now we suitably modify itg to obtain a sequence converging in L^{Q) to the linear function wo{y) = 
(Vu(xo),y). Set 

Ue{xo + gy) - u{xo) 


WeAv) '■= 


Ve,giy) ■= Veixo + QV)] 


then, letting first e —>■ 0 and then p —>• 0 we get {we,B,Ve,g) —>• {wq, 1) in L^{Q) x L^{Q). Moreover we 
have 


y-aixo) > liinlimsup / (anf (?/)|Vwg,g(i/)p + 

> lim limsup£’E(uig,g,Pe,g,Q), 
g^O £_).o 


P f {Ve,g{y) - 1)^ , I A ^ m 2 I I j 

Vl (-^-+ 7A»..<.(!/)I 11 iy 


where in the last inequality we have used that 1 /p > 1 for g small. 

By a standard diagonalization argument we can find two positive vanishing sequences (e/i), (p/t) such 
that —>■ (wojl) in L^{Q) x L^{Q) as ft. — >■ + 00 , and 

fJ-aixo) > liuiini 

n—¥-\-oo 

Let Q' CC Q] then, gathering Proposition 2.2(i) and Theorem 4.2, we get 


h—¥-\-oo 


Ma(^o) > lim inf J 


V, 


2 iv7,.. |2 I P (i'^sh,eh 1)^ 


£h,gh ^^h,eh 




2y/2 


£h 


+ c{Q,Q')el\V\eh^ 


Bh I 


dy 


> a|VM(xo)p/:"(Q'), 


and (5.2) follows letting Q' Q. 

We now prove (5.3). Let xo € S'^ be a Lebesgue point for g with respect to *-e., 


yj{xo) = lim 


ft((5e(2^o)) 


= lim 


ft(Q^(^o)) 


e^O-H"-l(Q^(a;o)nS'„) g^Q p"-i 


(5.6) 


where v := i'u{xq). The Besicovitch derivation Theorem ensures that (5.6) holds true for 'H” ^-a.e. 
xq & Su- By the definition of approximate discontinuity point, we can assume that in addition 


i 


lim 

e^O^(Q-(a:o))± 


i{x) — M^(a;o)| dx = 0 , 


(5.7) 


where {QY,ixo))^ := {x e (Qg(a;o)): ± (x - xq, v) > 0 }. 
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By following the same argument as before, we get 

If 1 

> 1 ™ limsup —^/ie(Q^;(a;o)) 

e" 1 JQ-ixn) P" ^ 


= lim lim sup 

£»”" 


-/ ( 


avl{x)\Vue{x)\'^ + - IL _|_£3|A^^(|a.)|2 


limlimsupp / ( auf (xq + p?/)|VMe(xo + py)| 

e^o £->0 JQ*- V 

/3 f {ve{xo + gy) - If 


2^2 

2 


dx 


+ 


2^2 


= lim lim sup 


^“^0 g_>o Jq'^ V 




■e^|Aue(xo + py)|‘ 


dy 


(I‘"“‘i'* - +(;) <'!' 

> lim limsup£■£/„(«£,(,, 
e-^-o e^O 

where we set Ue^g{y) := Ue(xo + gy), Ve,g{y) '■= Ve(xo + gy)- Notice that in view of (5.7), letting first 
e —>■ 0 and then p —>• 0 , we get {ue^g,Ve,g) —>■ (uo, 1 ) in L^iQ") x L^{Q'^) where 


uo(x) := 


u“''(xo) if (x — xo, v) > 0 , 
u~(xo) if (x — Xo, u) < 0 . 


Then, a diagonalization argument provides us with two positive vanishing sequences (sh), (gh) such that 

(5.8) 


cr/t := -)■ 0 as h -)> +oo, iue,^,g^,Ve^,g,f) -)■ (uo, 1 ) in L^{Q'^) x L^{Q'') as h ^ +oo, and 

A^j(3^o) ^ ^tyhf^eh,QhT'^Sh,QhTQ )• 

ri—¥-\-oo 

Set Uh '■= UehiBh^ '^h ■= '^eh,eh- 

Since is invariant under translations in u and under rotations in u and v, it is enough to bound 
from below the right hand side in (5.8) when u = e„ and Uh —>■ Ut := iX{a;„>o}; t ■= ■u“''(xo) — u“(xo). 

To this end let d > 0 and Q{S) := (—1/2 + d, 1/2 — d)"; then. 


/.‘j(a;o) > lim £ahiuh,Vh,Q) > lim £a^{uh,Vh,Q{6)), 

/i—^+oo h—^-\-(x> 


for every d > 0 small. We now show that 


lim £cr^ {uh, Vh, Q{S)) > P. 

h—¥-\-CX) 


(5.9) 


(5.10) 


The proof now follows the line of that of [24, Lemma 3.4] where the asymptotic behaviour of a variant of 
the Modica-Mortola functional is investigated. The idea is to estimate from below the functionals Ea^ in 
a way which allows us to reduce to the one-dimensional lower bound proved in Theorem 4.1. 

In order to not overburden notation we now drop the index h for the sequences of functions as in 
(5.9). We assume moreover that v G C°°{Q{5)) and we write 

^xV = v^z + AyV, y G (5'(d), z G (-1/2 -|- d, 1/2 - d), 
with Q'{5) := (—1/2 -I- d, 1/2 — d)"“^. Then, we have 


{u,v,Q{S))=f f fav'^ {ul + \Vyu\'^) 
J-^+sJq'{5) \ 


P ({v-lf 


2^2 




+ O'? \Vzz + A„u| 


dy dz. 


(5.11) 
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We are going to estimate the right-hand side of (5.11) from below with a functional that no longer contains 
partial derivatives with respect to y. Since the term involving Vyit is non-negative the only term that we 
have to estimate is the one containing For y G CI^(Q(S)), we write 


f f 

r I 


\vzz + Aywl ?7 dydz 

r rh-S 

,,2 J2 


v^^r] dydz + 


I-1+sJq'{S) 


[\Ayv\'^ + 2vzzV^ ^yv) dydz. (5-12) 


We first estimate the last term in (5.12) from below. We have 


r if 

J — -2+'5 Q' 

-f'~‘ f 

J-i+S JO' 


Vzz 'tf AyV dy dz 


i+SJQ'iS) 
-5 r 


Vz AyVz dydz — 2 


r! 

J-^+5JQ'{S) 


Vz rj Tjz AyV dy dz 


= / {iVyVzl’^r]'^+ 2vzV{^yVz,^yV)) dydz-2 

J-^+s Jq'{s) 

By Young’s inequality we get 

i-5 r f^-S f 

Vz-n {VyVz,'^yri) dydz < 


r‘ f 

J-^+S JQ’iS) 


Vz r] r]z AyV dy dz. 


2 

and 


r I 

J-^+S JQ'iS) 




r‘ f 


l-J+SJQ'iS) 


'^yVz\ V dy dz -I- 


Vz r] riz AyV dy dz < 


- r' f 

2 J-i+5 Jq'{S) 


Ay rip ? 7 ^ dydz + 2 


l-i+sJQ'iS) 


n 

J-^+sJQ'iS) 


vl \S/yy\^ dydz 


vl dy dz. 


Using in (5.12) the two bounds as above we obtain 


[ [ [lAyvl"^+2vzz Ayv) r]'^ dydz >- 

J-^+sJq'{s) 


[ vl{Ayl + 2\Wyr,\^) dydz. (5.13) 
Jq'(s) 


i-5 

-1+SJQ'{S) 


Thus, by (5.11)-(5.13), we get 
-i-i 


crt 


l-J+SJQ'iS) 

Hence, if we assume in addition that 

rh-s 


[ [ \vzz + Ayv\^y^dydz>al[ [ v^^y^ dy dz - c{y)al [ vldx. 

J-J+5JQ'{S) 7-1+5 Jq'(5) Jq{S) 


loo < 1, we can conclude that 




I_i+5 7q'(5) 


av \uz\ + 


/3 ({v- If 


2y/2 


(Jh 


+ <JhVzz] P7 dy dz 


2 V 2 


iv)<Jh [ ■ 

JQ{S) 


vt dx. 


which holds also true for v € lU^’^(H), by virtue of the density of C°°{Q) fl lU^’^(H) in lU^’^(r2). 
Appealing to Proposition 2.1 and Proposition 2.2(i) we have 


13 

2 V 2 ' 


[ iv)l 

Jq{s) 


dx < ca 


2 /3 


2\f2 yjQiS) 

< alc{Q,Q{S))£^y{u,v,Q)] 


+ cr^|V^i;p dx 
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therefore 

lim £ahiuh,Vh,Qi6)) 

h—>--\-oo 

- f + - — + crl{vh)l^'\'\r]^dydz, 

h-^+oo J\ 2V2 V JJ 

for every t] e C^{Q{6)). Hence if we choose S > 6 and 77 S C^{Q{6)) such that 77 = 1 in Q{S) CC Q{S), 
invoking Fatou’s Lemma and Theorem 4.1 we get 

dji^o) > lim £^^{uh,Vh,Q{S)) 

- / - f ~(^^hiuh)l +^ (— - — + al{vh)l,')') dz dy 

^-*■+ 0 ° Jq'( 5) J_i+5 V 2v2 \ (Jh JJ 

- f / (<^^hiuh)l +^ (— - — + al{vh)l^'\'\ dz dy 

jQ'iS) J-^+s \ 2 v 2 \ (Jh JJ 

> pn^-\Q'{5)). 

Then, (5.3) follows by letting 5 —>■ 0. 

If M G GSBV^{£l) the thesis follows by a standard truncation argument. In fact if := (it A M) V 
(—M), then € SBV^(fi) for all M € N. Then, appealing to the lower-semicontinuity of A4S and 
noticing that £e{-,v) (and hence £'{-,v)) decreases by truncation, we immediately get 

£'{u, 1) > liminf MS{u^, 1) > MS{u, 1). 

Af—^+00 

□ 


The following result holds true. 

Proposition 5.2. IFe have 

T-lim fe(it, 1) = AI5(it, 1) (5-14) 

e->0 

for every u G GSBV^(fl). 

Proof. The lower bound inequality is a consequence of Proposition 5.1 while the upper bound can be 
proved by taking the same recovery sequence as in Theorem 4.2. □ 


Now define 


£e{u,v) := 


£e{u,v) 

G IFo 

and vVu G T^(H; R"), 

+00 

otherwise in L^(H) x 


then we can prove the following T-convergence result. 

Theorem 5.3. The sequence {£e) T-converges, with respect to the {L^{£1) x L^{£1))-topology, to the 
functional A4S as in (2.5). 


Proof. The F-convergence result is a straightforward consequence of Proposition 5.14 once we notice that 
thanks to the boundary conditions satisfied by v, we can now invoke Proposition 2.2(ii) to get 

£e{u,v)> j ^ai!^|Vitp + —h c(H)e^|V^iip^^ dx, (5.15) 
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which together with (5.1) allows us to conclude that the domain of the F-limit is GSBV^{ft) x {v = 
1 a.e. in fl}. □ 


Remark 5.4. The C^-regularity of dil is only used to invoke Proposition 2.2(ii) in order to obtain the 
estimate from below (5.15). We notice however that for n = 2, which is the interesting case in numerical 
simulations, Proposition 2.2(ii) holds also true in bounded polygonal open sets (see e.g. [23, Theorem 
2.2.3]). 

5.1. Convergence of minimization problems. For every {u,v) S x VF^’^(fl) consider the 

functionals 

Ee{u,v) + 'qe / |Vupdx + 7 / \u — g\^ dx. (5.16) 

J Q. J Q, 

Appealing to Theorem 5.3 and to the fundamental property of P-convergence, also in this case we can 
prove a result on the convergence of associated minimization problems. 


Theorem 5.5. For every fixed e > 0 there exists a minimizing pair {ug^Vif) for the problem 
Mg := ini < £g{u,v) + r]e I \Vu\'^dx + j j \u — g\^ dx: {u,v) G x W^’^{El) 

L 'f O, J Q 

Moreover, up to subsequences, {ug,Vg) — > {u, 1) in x M{fl) where u is a solution to 



if n = 1 then u G SBV'^{Vt). Further, Mg -G M as £ —> 0. 


(5.17) 


Proof. The proof follows the line of that of Theorem 4.5 once we notice that the convergence (4.32) is 
now ensured by Proposition 2.2(ii). □ 


6. Numerical Results 

In this section we discuss some numerical results for the second-order approximation £g{u,v); i.e., we 
minimize (5.16) for reasons of practicability. First of all, the discretization of the Laplacian is more 
straightforward and leads to more compact schemes compared to the discretization of the Hessian 
involving mixed derivatives. As we have noticed in Remark 4.3 also the Hessian penalization can be 
rewritten into a Laplacian one under certain conditions, which are however not suitable for a numerical 
implementation since we need to enforce a solution in The latter means we have to implement 

simultaneous Dirichlet and Neumann boundary conditions, which is not feasible in finite difference or 
finite element discretizations without enforcing additional constraints. A second argument comes from 
the comparison with the first-order Ambrosio-Tortorelli functional, which already includes a Laplacian in 
the optimality condition, while the optimality condition for (5.16) changes only to a concatenation of two 
Laplacians. Hence, the modification of a code for the Ambrosio-Tortorelli functional to the second-order 
version (5.16) is straightforward and allows for a comparison of computational efficiency. 

We shall report on several computational experiments, starting with simple synthetic images that 
allow for a detailed study of fine properties such as the realization of the optimal profile already for a 
rather low number of pixels. Subsequently we investigate the behaviour on a set of natural and biomedical 
images, highlighting several differences of the second-order approach to the classical Ambrosio-Tortorelli 
functional. All experiments are carried out in two spatial dimensions, but we mention that extensions to 
volume data sets are obvious. Since a combination of the parameters a, /3 and 7 is redundant, we choose 
/? = 0.3 in all experiments. 
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6.1. Numerical Solution. In order to minimize the functional in (5.16) we follow the common strategy 
of iterative minimization. Hence, given an iterate we subsequently compute 

G argminfe(-«^, u) 

V 

G argminfe(M, -I-77e / \Vu\'^dx + j / \u — g\^dx. 

“ Jn in 

This yields a descent method for the overall functional, which is based on solving two quadratic minimiza¬ 
tion problems, respectively the corresponding linear optimality systems in each case. The same method 
is used to minimize the original Ambrosio-Tortorelli functional in (1.2). 

The linear equation to be solved for is in both cases 


-2aV • ((u'=+^)2Vu''+i) - 7j,Au'^+^ + = yg. 


The linear equation for is given by 

2a\Vu^\‘^v^+^ + 

\/2e 72 ^/2e 


in the case of the new second-order functional, respectively by 


2a|v777+^ -h -7+^ - /3eA7+^ = - 

e e 


in the case of the Ambrosio-Tortorelli functional. 

We discretize the functionals in (5.16) and (1.2) by standard finite differences on a rectangular grid, 
using one-sided (forward) differencing for the gradient and the adjoint one-sided differences for the diver¬ 
gence, hence the usual central differencing for the Laplacian. For the discretized linear systems we use a 
direct solver. Convergence diagnostics and stopping rules are based on the relative size of the change in 
the images; i.e., 


e 


k 


= max 


||7+1-7|U l|7+i-7||ool 

h'^+^lloo ’ |k'=+M|oo i 


6.2. Simple Test Examples. We start with a simple image showing a one-dimensional structure in the 
vertical direction. Figure 6 illustrates the results for the parameter settings a = 10 “^,7 = 10“^, e = 
3* 10“^, which yield a visually optimal result at the given image resolution. The segmentations obtained 
from the two models (see images of v in the middle) are not distinguishable by eye, a fine difference can 
be seen however when exploring the level sets where v is slightly larger than one in the second-order 
model. Note that due to a maximum principle respectively the monotone shape of the optimal profile 
in the Ambrosio-Tortorelli model the variable v is always less or equal one, while the optimal profile in 
the second-order model exceeds one, hence the corresponding level set shall provide further information 
about edge location. This is illustrated in the right-most plot in Figure 6 , from which one observes that 
V in the second-order model can provide an approximation of the edge set from both sides - an accurate 
reconstruction can be obtained as the midpoints between the local maxima. This behaviour is also present 
for larger values of £ as illustrated in the supplementary material. 

With the same set of parameters we also compute minimizers for images of an ellipse with large ratio 
between the main directions, illustrated in Figure 7 and two overlapping circles, illustrated in Figure 8 . 
Both images are perturbed by additive Gaussian noise to test also the effect of noise on the results. We 
again plot the resulting minimizers v for both models and the level set for the second-order models. The 
clean images u yield no visible differences and are shown in the supplementary material for completeness. 
Overall we observe analogous behaviour as for the one-dimensional example, a remarkable fact is that the 
level sets of v in the second-order model are able to provide a well separated segmentation of the overlap 
region. 
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Figure 6. One-dimensional structure (from left to right): Image g, resulting v in the 
Ambrosio-Tortorelli model, resulting v in the second-order model, binary plot of the level 
set {r; > 1.005} in the second-order model. 



Figure 7. Ellipse (from left to right): Image g, resulting v in the Ambrosio-Tortorelli 
model, resulting v in the second-order model, binary plot of the level set {v > 1.005} in 
the second-order model. 






Figure 8. Two circles (from left to right): Image g, resulting v in the Ambrosio- 
Tortorelli model, resulting v in the second-order model, binary plot of the level set 
{v > 1.005} in the second-order model. 

6.3. Natural and Biomedical Images. In the following we report on the results of the second-order 
model for few examples of natural images as well as different kinds of microscopy images. A general 
observation on all those images is that at the given resolution (rather small compared to the size of the 
structures) and the hence possible choices of er, the convergence to the optimal profile is far less pronounced 
than for the simple images above. Consequently the level set plot does not show an approximation of 
edges from both sides, for this reason we do not display the plots here. On the other hand we observe 
more interesting behaviour in the clean images u, which we also provide in the supplementary material. 
Let us mention that the F-convergence to the same minimizer does not mean that the original Ambrosio- 
Tortorelli functional and the second-order approach yield the same or very similar results on real images, 
which is due to many effects such as the given finite resolution of the image and choice of e, the details 
of convergence in e, as well as the level of convergence in the numerical minimization. 
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The algorithm was tested on various natural images such as the Kodak image test set (see also 
the supplementary material). The most pronounced difference between the original Ambrosio-Tortorelli 
model and the novel second-order version concerns structures at a small scale, which is however still larger 
than the typical scale of noise. This is well illustrated in a portrait photograph of a person containing 
freckles (see Figure 1, for parameters a = 3 * 10“^, 7 = 3 * 10“^,£ = 7 * 10“^ ). Such a behaviour is 
observed also for a wider range of parameters and seems clearly related to the stronger smoothing of 
the second-order model in higher frequencies. Another - at least visual - impression confirmed also in 
other results is that the contours being present in the results of both models appear smoother in the 
second-order model, which may be advantageous in many cases. 

We also report on an effect we obtain for rather large choice of e, i.e. rather far from convergence. 
This is illustrated in a phase-contrast microscopy image of a mitotic cell (cf. [21]). Such images are 
challenging for segmentation algorithms due to halo effects at the border of the cell, while the interior 
has similar grey value as the surrounding medium. Choosing e much larger than in the examples before 
(a = 3 * 10“^, 7 = 3 * 10“^, £ = 5 * 10“^) we obtain a contour v that actually fills the whole interior. This 
effect is clearly benefitial for the second-order model, whose result allows a simple tracking of the cell. 



Figure 9. Mitotic cell, from [21] (from left to right): Image g, resulting v in the 
Ambrosio-Tortorelli model, resulting v in the second-order model. 

We finally comment on the convergence behaviour of alternating minimization algorithms routinely 
used for approximating minimizers of the Ambrosio-Tortorelli functional. The observation made in the 
majority of our numerical experiments is that the number of iterations needed for fixed accuracy in the 
second-order model is at least comparable to those for the standard first-order version, in many cases 
the number of iterations is significantly reduced for the second-order model in particular for real images. 
In some parameter cases a visual comparison of results indicates that one obtains a global minimizer 
for the second-order model, while the iteration for the first-order model is stuck in a suboptimal local 
minimum, which is however difficult to verify. We refer to the supplementary material to a collection of 
computational investigations of convergence. 
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Appendix A. Supplementary material 

A.l. Convergence Behaviour of Alternating Minimization Algorithms. The convergence indi¬ 
cator is plotted vs. the number of iterations in several examples in Figure 10 and 11. Figure 12 
illustrates two cases of parameters where alternating minimization on the Ambrosio-Tortorelli functional 
did not converge. 



Figure 10 . Convergence history of vs. number of iterations k. One-dimensional 
example (left, a = 10 “^, 7 = 10 “^, e = 9 * 10 “^), ellipse (middle, a = 10 “^, 7 = 10 “^, e = 
3 * 10 “^), two circles (right, a = 10 “^, 7 = 10 “^, e = 3 * 10 “^). 



Figure 11 . Convergence history of vs. number of iterations k. Kodak image nr2 
(left, a = 10 “^, 7 = 10 “^, e = 3 + 10 “^), nr 7 (middle, a = 10 “^, 7 = 10 “^, e = 7* 10 “^), 
nr 23 (right, a = 10 “^, 7 = 10 “^, e = 7* 10 “^). 


A.2. One-dimensional Structure. Figure 13 displays further results for the one-dimensional structure, 
parameter a = 10“^, 7 = 10“^, c = 9 * 10“^. 

A.3. Results on Kodak Images. We display some examples of results on the Kodak images 2, 7, and 
23, displayed in Figure 14. Figures 15, 16, 17 display the resulting segmentation v for the Kodak image 
nr 2 with a = 10“^ and different values of 7 and e. Figures 18 and 19 display the results for the Kodak 
image nr 7 with e = 7 * 10“^ and different values of a and 7 . Figure 20 displays the resulting v in Kodak 
image nr 23 for a = 10“^, 7 = 10“^, e = 7 * 10“^. 

A.4. Reconstructed Images. Figures 21, 22, 23, and 24 display the resulting u in the different models 
with the parameter settings in the paper. Figures 25, 26, 27 display results for the Kodak image nr 2 
with a = 10“^ and different values of 7 and e. Figures 28 and 29 display the results for the Kodak image 
nr 7 with e = 7 + 10“^ and different values of a and 7 . Figure 30 displays the resulting u in Kodak image 
nr 23 for both models with parameters a = 10 “^,7 = 10“^,e = 7* 10“^. 
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Figure 12. Convergence history of e^' vs. number of iterations k. Kodak image nr 7 
(left, a = 10“^, 7 = 7 * 10“"^, e = 7 * 10“^), Kodak image nr 23 (right, a = 10“^, 7 = 
7* 10-4,e = 7* 10-2). 



Figure 13. One-dimensional structure (from left to right): Image g, resulting v in the 
Ambrosio-Tortorelli model, resulting v in the second-order model, binary plot of the level 
set {t! > 1.005} in the second-order model, e = 9 * IQ-^. 



Figure 14. Kodak image nr 2 (left), nr 7 (middle) and nr 23 (right). 
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Figure 15. Kodak image 2: resulting v in the Ambrosio-Tortorelli model (left) and in 
the second order model (right), 7 = 10“^, c = 3 * 10“^. 



Figure 16. Kodak image 2: resulting v in the Ambrosio-Tortorelli model (left) and in 
the second order model (right), 7 = 7 + 10“^, c = 6 * 10“^. 
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Figure 17. Kodak image 2: resulting u in the Ambrosio-Tortorelli model (left) and in 
the second order model (right), 7 = 7 + 10“'*, e = &* 10“^. 



Figure 18. Kodak image 7: Resulting v in the Ambrosio-Tortorelli model (left) and in 
the second order model (right), both with a = 10 “^, 7 = 7 + 10“^. 
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Figure 19. Kodak image 7: Resulting v in the Ambrosio-Tortorelli model (left) and in 
the second order model (right), both with a = 1 * 10 “^, 7 = 10 “^. 



Figure 20. Kodak image 23: resulting v in the Ambrosio-Tortorelli model (left) and in 
the second order model (right). 
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Figure 21. Ellipse: resulting u in the Ambrosio-Tortorelli model (left) and in the second 
order model (right). 



Figure 22. Two circles: resulting u in the Ambrosio-Tortorelli model (left) and in the 
second order model (right). 
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Figure 23. Sisse image: resulting u in the Ambrosio-Tortorelli model (left) and in the 
second order model (right). 



Figure 24. Mitosis image: resulting u in the Ambrosio-Tortorelli model (left) and in 
the second order model (right). 
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Figure 25. Kodak image 2: resulting u in the Ambrosio-Tortorelli model (left) and in 
the second order model (right), 7 = 10“^, e = 3 * 10“^. 



Figure 26. Kodak image 2: resulting u in the Ambrosio-Tortorelli model (left) and in 
the second order model (right), 7 = 7 + 10“^, e = & * 10“^. 












SECOND-ORDER AMBROSIO-TORTORELLI FUNCTIONAL 


39 



Figure 27. Kodak image 2: resulting u in the Ambrosio-Tortorelli model (left) and in 
the second order model (right), 7 = 7 * 10“'*, e = & * 10“^. 



Figure 28. Kodak image 7: Resulting u in the Ambrosio-Tortorelli model (left) and in 
the second order model (right), both with a = 10“^, 7 = 7* 10“^. 
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Figure 29. Kodak image 7: Resulting u in the Ambrosio-Tortorelli model (left) and in 
the second order model (right), both with a = 7 * 10“^,7 = 7* 10“^. 



Figure 30. Kodak image 23: resulting u in the Ambrosio-Tortorelli model (left) and in 
the second order model (right). 
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